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MULTIPLICITY OF SOLUTIONS FOR LINEAR PARTIAL 
DIFFERENTIAL EQUATIONS USING (GENERALIZED) ENERGY 

OPERATORS 

JEAN-PHILIPPE MONTILLET 


Abstract. Families of energy operators and generalized energy operators 
have recently been introduced in the definition of the solutions of linear Par¬ 
tial Differential Equations (PDEs) with a particular application to the wave 
equation ESI- To do so, the author has introduced the notion of energy spaces 
included in the Schwartz space S~(R). In this model, the key is to look at 
which ones of these subspaces are reduced to {0} with the help of energy opera¬ 
tors (and generalized energy operators). It leads to define additional solutions 
for a nominated PDE. Beyond that, this work intends to develop the concept 
of multiplicity of solutions for a linear PDE through the study of these energy 
spaces (i.e. emptiness). The main concept is that the PDE is viewed as a gen¬ 
erator of solutions rather than the classical way of solving the given equation 
with a known form of the solutions together with boundary conditions. The 
theory is applied to the wave equation with the special case of the evanescent 
waves. The work ends with a discussion on another concept, the duplication 
of solutions and some applications in a closed cavity. 


1. Introduction 

The energy operator was initially called the Teager-Kaiser energy operator m 
and the family of Teager-Kaiser energy operators in [12j . It was first introduced in 
signal processing to detect transient signals [6] and filtering modulated signals [3]. 
After three decades of research, it has shown multiple applications in various areas 
(i.e. speech analysis[B], transient signal detection[5], image processing [2, optic 
[17], localization [lH|). This energy operator is defined as in [T3|, [14], [15] and 
through this work. [13] defined the conjugate operator in order to rewrite the 
wave equation with these two operators. Furthermore, El defined the family of 
energy operator ('k/'j^gz and (ik^)fcgz in order to decompose the successive deriva¬ 
tives of a finite energy function /" (n in Z+ — {0,1}) in the Schwartz space. [15] 
introduced the concept of generalized energy operators. In the same work, it was 
shown that there is a possible application of the energy operators and generalized 
energy operators to define new sets of solutions for linear PDEs. 

This work is the sequel of [15]. It intends to define the concept of Energy Spaces, 
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which are sets (in the Schwartz space) of soiutions for a nominated iinear PDE. 
Theorem 2 states the mechanism for the functions /” and 9^/" to be soiutions of 
a nominated iinear PDE. Then, the work is extended to the case when the energy 
operator famiiies (appiied to /, '!'+(/), 'I'p (/)) and generaiized energy operator 
famiiies can aiso be soiution of the same equation via the Corollary 1. An overview 
of the concept in the particular case of the wave equation is to consider these addi¬ 
tional solutions as waves generated by a given PDE (i.e. the d’ Alembert operator 
for the wave equation □(/) = 0 [3]) with lower energy. This approach differs with 
the traditional way of solving a PDE using boundary conditions with a known form 
of the solutions. When using the definition of energy space (Definition 3, [IS]), the 
formulation of finding each solution (e.g., /", 9^/", ('I'i^(/))", ('^'f (/))") can 

be reformulated in looking into which energy space is not reduced to {0}. Section 
6 of [15] is the proof of concept using the wave equation and the evanescent waves 
as special type of solutions. The last section concludes with a discussion on the 
potential of using this model with non-linear PDEs. 


2. Preliminaries 

2.1. Energy Operators and Generalized Energy Operators in S“(IR). Through¬ 
out this work, /” for any n in Z+ — {0} is supposed to be a smooth real-valued and 
finite energy function, and in the Schwartz space S“(M) defined as: 


S-(K) = {/eC°°(M), sRp,<o|t'=||9^/(t)| <oo, VA:GZ+, Vj G Z+} 

( 2 . 1 ) 

Sometime /" can also be analytic if its development in Taylor-Series is necessary 
(e.g, the application to closed cavity in Section 4.31. The choice of /"■ (for any n 
in Z+ — {0,1}) in the Schwartz space S“(M) is based on the development in [14] . 
Section 2, because we are dealing with multiple integrals or derivatives of /" when 
applying the energy operators (4'^)fcgz+, (4'^)fcgz+ and later on the generalized 
energy operators. 

In the following, let us call the set J^(S“(M), S~(IR)) all functions/operators defined 
such as 7 : S“(K) —?■ S“(K). Let us recall some definitions and important results 
given in [14] and m- 


Section 2 in m and Section 4 in (TS] defined the energy operators 4'^ {k in 
Z) and the generalized energy operators and (p in Z+). [15] defined: 



= dfdt 


= 


= 


k gZ 
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ipj, is the operator conjugate of Furthermore, m defined the generalized 
energy operators and : 




= atik+(.)5f-ivi-+(.) + vi/+(/)afvi'+(.) 


+ 

+ 

= a*[[.]°]+ari[[.]0]+ + [[.]°]+a^=[[.]°]+ 
= IIVU 


’ -Ife ]fe 



’ 'l/c ]fe 





( 2 . 2 ) 

By iterating the bracket [■], m defined the generalized operator and the 

conjugate [[■]^]'^ with p in Z+. Note that = 0 for all p in Z+. 

In addition, the derivative chain rule property and bilinearity of the energy oper¬ 
ators and generalized operators are shown respectively in El, Section 2 and El, 
Proposition 3. 

Definition 1 [T3]: For all / in S~(IR), for all v G Z+ — {0}, for all n G Z+ and 
n > 1, the family of operators (^k)kez (with (^k)kez Q -7^(S~(I^)) S“(ffi))) decom¬ 
poses 9)'/" in K, if it exists (-^j)jez+u{o} ^ Q K, and it exists (aj) 

and I in Z+ U {0} (with I < v) 

such as du- = y:~ 1 T.kL-N, Ck^k{d‘^>‘f). 

In addition, one has to define s“ (M) as: 

s-(K) = {/ e S-(K)|/ i (Ufcezifer(vk+)) U (Ufcez-{i}^fer(il'-))} (2.3) 

Ker{^^) and Ker{'^~^) are the kernels of the operators and (k in Z) (see 
m, Properties 1 and 2). One can also underline that s (K) C S (K) Following 
Definition I, the uniqueness of the decomposition in s“(IR) with the families of 
differential operators can be stated as: 

Definition 2 [T3]: For all / in s“(]R), for all v G Z+ — {0}, for all n G Z+ and 
n > f, the families of operators ('i'^)kez and ('i'j:)kez ((^'fc)feGZ and ('l'j:)kez Q 
J^(s“(IR), S~(IR))) decompose uniquely /" in K, if for any family of operators 
{Sk)k£Z C J^(S“(IR), S~(IR)) decomposing /" in M, there exists a unique couple 
in such as: 

Sk{f)=PiKif)+P2'i^~kif)^ VfceZ (2.4) 


Two important results shown in [14] (Lemma and Theorem) are: 

Lemma 0: For / in S“(]R), the family of DEO {k = {0, ±1, ±2,...}) decom¬ 
poses the successive derivatives of the n-th power of / for n G Z"*" and n > f. 

Theorem 0: For / in s“ (K), the families of DEO and (fc = {0, ±1, ±2, ...}) 
decompose uniquely the successive derivatives of the n-th power of / for n G Z+ 
and n > f. 
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The Lemma 0 and Theorem 0 were then extended in m to the family of gener¬ 
alized operator with : 


Lemma 1: For / in (R), p in Z+, the families of generalized energy operators 
(fc = |0,±1,±2, ...|) decompose the successive derivatives of the n-th power 
of [[fF~^]t for n e Z+ and n > 1. 

Theorem 1: For / in s“(IR), for p in Z+, the families of generalized operators 
and {k = {0, ±1, ±2,...}) decompose uniquely the successive deriva¬ 
tives of the n-th power of for n G Z+ and n > 1. 


S~(M) and s” (M) (p in Z+) are energy spaces in S~(IR) defined in the next section. 
Note that as underlined in [T3] (Section 3, p.74) and [TS] (Section 4), one can ex¬ 
tend the Theorem 0, Theorem 1, Lemma 0 and Lemma 1 for /” with n in Z. 
The appendix A recalls the discussion in [15] (Section 4). Here, n is restricted to 
Z+ — {0,1} in order to easy the whole mathematical development. 

Some time, the finite energy function /" (n in Z+ — {0,1}) can also be considered 
analytic. In other words, there are (p,g) {p > q) in such as /" can be developed 
in Taylor Series m- 

OO / \ /c 

rip) = + 

k^l 

(2.5) 


Proposition 1, m states: 

Proposition 1; If for any n G z+, r G s-(M) is analytic and finite energy; for 
any (p,q) G (with r in [q,p\) and £(/") in S~(R) is analytic, where 

£{r{r)) = f\r{t)ydt < oo ( 2 . 6 ) 

J q 

then 

£{np)) = £{nq)) + Y.df{nq)y^^£j^<^ 

k^O 

(2.7) 


is a convergent series. 

This property is specifically used in the applications of this work described in Sec¬ 
tion |4| 


2.2. Energy Spaces. The energy spaces were introduced in m, Definition 3 and 
m, equation (24). The definition reads: 

Definition 3: The energy space Ep, with p in Z+, is equal to Ep = UiGZ+u{o} 


with M* C S (R) for i in Z’*' defined 

M* = {5GS-(R)| g = 5*([[/f]+)",[[/f]+GS-(R), n G Z+- {0,1}} (2.8) 
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If g is a general solution of some linear PDEs, then /" can be identified as a special 
form of the solution (conditionally to its existence). One can further define the 
subspace S“(IR) C S~(K) for p in Z+ (see [T5] . equation (25)): 

S;(K) = {Ep= U MV{0}} (2.9) 

iGZ+U{0} 

The energy space Ep is said associated with £( [[•]^~'’^] )*"). Note that Ep is not 

empty, because the assumption is that ([[/]^]^)" is a solution of a given linear 
PDE throughout this work. Thus, the energy space cannot be defined without a 
nominated PDE. Now, one can define the subset s” (M) defined as: 

Sp-(K) = {/ e S;(K)|/ ^ U i'Jkez-{i}Keri[[fr]k))} (2-10) 

The subset s”(M) is also defined such as Ep ^ {0}. Thus, one can see that s” (K) C 
Sp (K). Note that S()'(]R) C s“(IR), because 8“ (K) C S“(IR) (pinZ+). Furthermore, 
one can also define the special case: 

N* = {geS-(K)| g = a*(/)",rGS-(M), n G Z+- {0,1}} (2.11) 

Following Section 4 in [15], E = UiGZ+u{o} ^ ® £ S“(IR). Finally, the 

energy space E is said associated with £( [[.]°] (*"). 


3. Theorems of Multiplicity of the Solutions for a given PDE 


3.1. Multiplicty of the Solutions. To recall the first section, a possible appli¬ 
cation of the theory of the energy operator is to look at the solutions for various 
values of n and i instead of solving the equation for specific values (e.g., boundary 
conditions). According to Equation (2.8) and Definition 3, it is equivalent to find 
which subspace M* {i G Z+) of Ep {p in Z+) is reduced to {0} (or respectively N® 
of E). Thus using Theorem 0 and Theorem 1, one can use the energy operator 
in order to find the subspaces N® or M® reduced to {0}. For example, let us define 
for i in Z’*' 


L® = {g G S- 


g = d\f = d^r\^U!) + 'f'ra)),/ G s- 


(3.1) 


with Equation (2.8), L® C N®. If |4')^(/)| = 0 then L® = {0}. Using Definition 1 
and Theorem 0, one can write if it exists ii in Z+ such as |9(®4'(j^(/)| = 0, then 
L®® = {0}. Subsequently for all L®^ = {0}. 


3.2. Statement of the Theorem on Multiplicity of the Solutions for a 
linear PDE. If / is a solution of a linear PDE, let us call this statement A^/ = 0. 
One can summarize the philosophy of multiple solutions generated by a linear PDE 


Theorem 2 : / in s (M). Then, dlf^ (i in Z+, n in Z+ — {0,1}) is solution for (t, 
r) in [a, 6]^ (a < b, (a, b) in M^) if it is assumed 

1. (general condition to be a solution) A^9j/®®(r) = 0 

2. (Solutions in S“(K) ), dlf^ is a finite energy function such as £{dlf^){T) < 
oo 

3. (3 2) it exists nii in K, for i in Z+ such as rrii = sup{vtigz+-{o i} TG[a &]} 

(£(9®r)(r)) 
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4. (Superposition of solutions and energy conservation ) F{t) = Efe6Z+ 
then £{F{t)) < oo 

5. (5 ^ 4) 3 ii in Z+ such as V z > Zi, N* = {0} 

From the statement (5.), one can then dehne the energy space E = lJi6[o p-i] 
Now, if we want to use the decomposition of the solutions with the energy 

operators (e.g, Theorem 0) in order to find Zi such as is reduced to {0}, then 
the definition of the subspace N* reads 

N* = {geS-iR)\ g = = + 

fes-{R), nez+-{0,l}, a„eK} (3.2) 

That is why we define Theorem 2 with / in s“(IR). Furthermore, Theorem 2 
can be extended to the generalized operators as solutions of the linear PDF. With 
p in and the energy space definition (e.g, Definition 3), the Corollary 1 of 
Theorem 2 reads 

Corollary 1: For p in Z+, / in s“(M), 9^([[/J^’])'')" (z in Z+, n in — {0,1}) is 
solution for (t, r) in [a, 6]^ (a < b, (a, b) in if 

1. (general condition to be a solution) for z in TL^ and rz in Z+ - {0,1}, 

= 0 

2. (Solutions in S (M) ) for i in Z+ and n in Z+ — (0,1}, i9(([[/]p]{’)" is a 
finite energy function such as £(9 {([[/]p]{')”)(t) < oo 

3. (3 2) it exists rrii in M, for i in Z+ such as nii = szzp{vnez+-{o i}} 

imi[[mtr)ir)) 

4. (superposition of solutions and energy conservation ) 
nr) = Efca+5'=([[/F]+)”)(r)<oo 

5. (5 ^ 4) 3 zi in Z+ such as V z > zi, M* = {0} 

To recall the remark at the end of the statement of Theorem 2, one can justify that 
/ in s“ (M) in order to use the decomposition with the generalized energy operators 
(e.g. Theorem 1). Thus, the definition of the subspace M® reads 

M® = {geS-{R)\ g = dl{[[fr]^)^ = 

an{di-\[[fr]tT-\[[fr^]t + [[/f+f:). 

/ € Sp (M), n € Z’*' — jo, 1}, an € K} 

(3.3) 


Theorem 2 and Corollary 1 could then be applied to a given linear PDF by re¬ 
placing the first statement (A^(.) = 0). The next section highlights the assumptions 
behind each assertion in Theorem 2 and Corollary 1. 

3.3. Underlining Hypothesis of the Theorems and Proofs. This section dis¬ 
cusses each statement in Theorem 2 and Corollary 1. We are not attempting 
to give a formal proof, because we conjecture the existence of solutions for a given 
PDF such as the energy spaces S' (M) and s“ (M) (p in Z+) are not reduced to {0}. 
Keeping that in mind, one can give the following explanations to understand the 
previous section. 
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3.3.1. Theorem 2. 

Proof. 1. This is the definition of a solution for a nominated PDE. 

2. All the solutions considered in this work are finite energy functions in order 
to apply the decomposition using the energy operators stated in Theorem 
0 and Theorem 1. Note that it does not rule out that for some solutions 
(particular values of i and n), one can have {£{dlf'^){T)) ~ oo. In this case, 
the solutions cannot be accepted. 

3. The existence of the upper bound of £(9(/")(r) ( r in [a, 6], Vn S — 
{0,1}) comes from the definition. However, let us define the subspaces Hi 
(i in Z+), Hi CM. such as 

H, = {ze K |0 = £{dir){T) < oo, r G S(M-), n G Z+ - {0},t G [a, b]} (3.4) 


One can then define rui in M [5] such as 

m, = sup{VTiGZ+-{0,i}.rG[a.b]}(^('9t/")('r)) (3.5) 

Furthermore, we can also define M such as 

M = maa;vigz+-{o}(wi) (3.6) 

Then, 

M = SMp{VnGZ+-{0,l},TG[p,il],ViGZ+}(^(5t/”)('^)) (3-7) 

4. This statement is to guarantee that there is a finite sum of energy. In other 
words, there is no infinite number of solutions. Thus with the development 
in statement (3.), one can use the Minkowski inequality (e.g, [7], Theorem 
202) for T in [a, b\ 

= f \Y,dir{t)\^dt 

iGZ+ 

(f(F(r)))°-^ < ^ 

iGZ+ 

(3.8) 

iGZ+ 


Thus, (4.) stands if and only if X]iGZ+ < oo. As for all j G Z+, mi is 
in it then exists io in TP such as Wi > io, then rrii = 0. 

5. Following the discussion in Section 2.2 (e.g. Definition 3), one can see that 
if it exists A such as for all f > A, N® = {0} then £{dlf‘^) — 0. Or with 
the above development in (4.), the sum of coefficients m° ® becomes finite: 


{£{F{T))f^ 

zGZ+ 


{£{F{T))f^ 

Zl-1 

< ' 

2=0 

(3.9) 


□ 
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3.3.2. Corollary 1. 


Proof. The structure of the Corollary 1 is very similar to Theorem 2. In order 
to avoid repetitions, some explanations are shortened. For p in Z+, 

1. This is the definition of a solution for a nominated PDF. 

2. All the solutions considered in this work are finite energy functions in order 
to apply the decomposition using the energy operators stated in Theorem 
0 and Theorem 1. 

3. Similarly as in the previous discussion, the existence of the upper bound of 

^ ill )(''') ( ''' ill [®> G ~ {0> 1}) comes from the definition. 

Similarly to the statement (3.) of Theorem 2, one can define Hf ii in 
Z+), iJf C K such as 

Hi = {z G R\z = £idmr]tr)ir) < oo,r G S(K-),n G Z+ - {0},r G [a,6]} 

(3.10) 


Subsequently, with rrii^p in K 

= SMP{VnGZ+-{0,l}.TG[a.b]}('S’(l9t([[/]n^)")(T)) (3.11) 

and defining M such as 


Mp = TOaxviez+-{o}("iGp) 


Then, 


(3.12) 


Mp = SUP{VnGZ+-{0,l}.re[p.9].VieZ+}('^^(5t([[/]^]f)”)M) (3-13) 


Following statement [3.] and the previous discussion, one can use the 
Minkowski inequality (e.g, [7], Theorem 202) for r in [a,b] 


£{F{t)) = 

(£(F(t)))°-' < 



IE mfimtri^dt 

iGZ+ 


E (r mifimtrM°-^ 

iGZ+ da 


iGZ+ 


(3.14) 


As previously underlined, (4.) stands if and only if < c>o. In 

other words, it then exists io in Z"*" such as Vi > io, then rrii^p = 0. 

5. Following previous discussions for Theorem 2 and the above development, 
the sum of coefhcients ® becomes finite: 

{£{F{T))f^ < J2<P 

iGZ+ 

(£(F(t)))°-^ < E<P (3-15) 

i^O 

□ 
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3.4. Discussion on non-linear PDEs. The current work started with the pre¬ 
liminary study in m where the energy operators showed how to decompose the 
Helmholtz equation. The Helmholtz equation was chosen as an easy application, 
because of its linearity and the form of the solutions. In [15j , the preliminary results 
were extended. However with the additional development in this current work, the 
introduction of the energy space with the definition based on the solution of a nom¬ 
inated PDE does not justify anymore the assumption that the PDE must be linear. 
In addition, the assumption on the linearity of the nominated PDE in Theorem 
2 and Corollary 1 does not play any role. Thus, one should also be able to apply 
this model to non-linear PDEs. 

4. Application to linear PDEs 

4.1. Comments on functions of two variables solutions of linear and Non¬ 
linear PDEs. In this section and the remainder of this work, the finite energy 
functions of one variable described in Section |2.1[ are now functions of two vari¬ 
ables referring to the space dimension (r) and time (t). One has to add in the 
notation of the operators the symbol t or r to indicate which variable the deriva¬ 
tives refer to. For example, the operators and (fc in Z, p in Z+) 

refer to their derivatives in space, whereas and [[.]^]^’* {k in Z, p in Z+) 

refer to their derivatives in time. Following the discussion in [TS] (Section 6), one 
can then define the Schwartz space S~(M^) for function of two variables such as: 

S-(K2) = {/ e c-(R), V(ro,to) e K+l sup,^o\t'^\\dif{ro,t)\ < oo, 
and sup^^f)\r'^\\d^^f{r,to)\ < oo, Vfc e Z"*", Vj e Z+} 

Following this definition, the extension of the subspace s“(M^) C S“(M^) is: 

s-(r 2) = {/es-(K2)|vfcez, vE-+’‘(/)^{0} 

VA:eZ-{l}, 4/-’‘(/)^{0}}U 

{/eS-(K2)|vfcez, vi/+’’'(/) ^ {0} 

VA:eZ-{l}, ^ {0}} 

In |15j . it was emphasized that Definition 0, Definition 1, Lemma 0, Lemma 1, 
Theorem 0 and Theorem 1 can easily be extended to the function of two variables 
using the above definitions of s“(R^) and S“(M^). We will not state formally all 
the work previously done for the case of the functions of two variables in S“(K^). 
It is only a matter of replacing the variables from time to space. 

According to the previous section, one can state in the case of a function of two 
variables that it exists (a^, an 2 ) in such as for / G s“(IR^), n G Z+ — {0,1}: 

I dir= anif5ri(r-2(vi/‘’+(/) + vi/i’-(/)))), 

[ d;r= + (/)))) 

The definition of the energy space E (e.g, Definition 3) is defined in S“(IR^) 

Nj = {gGS-iR^), ytGR+, roGmg{ro,t) = diriro,t), 
fGs~{M.^), n G Z+— {0,1}} 

N* = {gGS-{R^),yrGR, to G R+\g{r,to) = dirir,to), 
/Gs“(M^), n G Z+— {0,1}} 


(4.2) 
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Thus, E = {U^ez+u{o} ^ {0}}U{Uiez+u{o} ^2 + {0}}, ECS (M^). Similarly, 
one can further extend the definition of energy space with the generalized energy 
operators in 8“ (M^) (p in Z+) with defining and M^. 

Furthermore, it is worth underlining that if / is a solution of a given linear PDE, 
then /"(ro,<) is in and /"(r, <o) is in N^. Choosing a wise choice for the 
definition of / (i.e. amplitude not equal to 0 ) leads to N? and N 2 not equal to { 0 }. 
One can then conclude that the energy space E ^ {0} and E 7 ^ {0} . 

4.2. Evanescent waves and the wave equation. The evanescent waves were 
already chosen in the numerical example in Section 6 of na. Here, this type of 
solutions of the Helmholtz equation is used as an application of the multiplicity 
theorem stated in Theorem 2. From m or [I], the Helmholtz equation can be 
formulated ; 


9fgir,t) 
or □^(r, t) = 0 , 
t G [0,T], r G [ri,r 2 ], 
to G [0,T], ro G [ri,r 2 ] 


-j2d?gir,t) = 0 , 

(ri,r2,T) G ri < r2 


(4.3) 


c is the speed of light. Note that the values of t and r are restricted to some interval, 
because it is conventional to solve the equation for a restricted time interval in M"*' 
and a specific region in space. Furthermore, it is well-known that the general 
solution g{r, t) of this equation is a sum of two waves traveling in opposite direction 
such as g{r, t) = fi{t — r/c) + fi{t + rjc) (e.g., [E]). With the development in the 
previous sections and taking the wave traveling along the increasing positive axis 
(if the other solutions is considered to travel along the increasing negative axis), we 
are interested in the solutions of the kind g{r, t) = i9)/"(r, t) (n in Z+ — {0,1}, p in 
Z+). 

As underlined in m. solutions in S (M^) of equation ( |4.3[ ) are finite energy 
functions such as the ones decaying for large values of r and t. This is a very 
limiting condition. For example, planar waves are not included. 

Now, let us apply the multiplicity theorem (e.g. Theorem 2) to the evanescent 
waves [ 1 ] : 


= 

t G [0,T], r G [r-i,r 2 ]. 


i?eaZ{Aexp (fc 2 r) exp (j(wt — fcir))}, 
(ri,r2) G ri < r2 


(4.4) 


ki and ^2 are the wave numbers, uj is the angular frequency and A is the amplitude 
of this wave m- Assuming w and (fci, ^ 2 ) known, one can add some boundary 
conditions in order to estimate fci, k 2 and A. However, our interest is just the 
general form assuming that all the parameters are known. 

The solutions in Nj and N 2 can be stated 


diriroA) = 
dir{r,to) = 

t G [0,T], r G [ri,r 2 ], 
n G Z+ - {0,1}, f G Z+ - {0} 


(jwn)V”(ro,t)}, 

(n(/c 2 - jfci))V”(r, to)}, 
(ri,r2,T) G ri < r2 

(4.5) 


to G [0,T], ro G [ri,r 2 ], 


It is well-known that a traveling wave solution of equation ( |4.3[ ) should satisfy the 
dispersion relationship between ki, k 2 and w [16]. Note that t is in G Z+ — {0}, 
because we are doing an application where the energy operators should help to 
determine which energy space is not empty. For i equal 0, the energy operators are 
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Table 1. Parameters for the scenarios 1 and 2 


Parameters 

Scenario 1 

Scenario 2 

ki {m~^ ) 


0.001 

h ) 



A (to) 

1 

1 

uj ( sec“^) 

0.03 

3000 

T ( sec) 

100 

100 

[ri,r2] ( to) 

[2,140] 

[2,140] 


not involved. This space can be seen as the form of desired solutions (with n in 
Z+-{0,1}). 

If we replace the general solution g in equation (4.3) with solutions in and N 2 
(see Section 4.1), one can study the dispersion relationship for different solutions 
of the form dlP{r,t) or dlP{r,t) : 


□5)r(r,t)=0 
Real{ — -^-(n.(/c 2 — jfci))^} = 0 

Real{^^^ - {ik2 - jh))'^} = 0 

□a;r(r,t) = o 
Real{ — -^-(n(fe — jki))^} = 0 

Real{^^^ -((/c2 - jfci))^} = 0 (4.6) 

c 

It means that the dispersion relationship for this type of solutions is a function 
neither of the degree of the derivatives i nor the power n for the special case of the 
evanescent waves. Furthermore, one can also calculate the closed-form expression 
9jT)'"’*(r, t) and i9).4')^’’'(r,t) such as 

= Real{{j2iu})f{ro,t)} 

= Real{{2i{k2- jki))f{r,to)} 
tG[0,r], rG[r’i,r2], (ri, r2, T) G ri < r2 

toG[0,r], roG[ri,r 2 ], nGZ+-{0,1}, i G Z+- {0} (4.7) 


These expressions can help to determine which subspaces and N 2 can be re¬ 
duced to {0} according to the model explained in Section 3 (e.g. Theorem 1). 
Now, let us estimate numerically ii such as = 0 and 12 such as 

|g 22 \j>+T(j)| _ Q Jqj. fjxed values of A, ki, k 2 , w (in K) and a given region in 
space and interval in time. Note that it is mathematically more accurate to state 
|gn^j>+.t(j)| ^ Q |9*2T)''’’'(/)| ^ 0 rather than using the symbol ’=’, because we 
are estimating numerically the different solutions. Thus, one can decide implicitly 
that it exists e in K such as if |9j^4')'’’*(/)| < e , then 19)^4')'’’*(/)| = 0 (reciprocally 
|5;i^^’’(/)| < e , then |9;iT+’’'(/)| = 0). 

The parameters displayed in Table 1 are for two scenarios. The parameters in 
Scenario 1 are used in the simulation to estimate ii and reciprocally Scenario 2 for 
i 2 - As previously mentioned in the definition of the energy spaces (e.g. Definition 
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Figure 1. Results from the estimation of the energy of 
and t)). The Figures correspond to 

Scenario 1. 






3), the energy function is said "associate" with the energy space. In other words, 
one can estimate the energy of t) and to) for the various 

values of i in Z+ — {0}. Figure 1 and 2 displays the results. Note that the simula¬ 
tions are performed using Matlab (v9.0 R2016a - mathworks.com). 

The top figures in Figure 1 and 2 correspond to the numerical estimation of 
the energy of the operators {f){rQ,t) (also called t))) and f 

9).il']'’’’'(/)(r,to) (also called 9*[[/]°])'’’'^(r,to))). Thus, the results shows that for 
t > 2, the values of £(i9j4')'’’*(/)(ro, t)) and £’(i9).4')^’’^(/)(r, to)) are smaller than 
IQ-io^ Thus, we can conclude that ti = ^2 = 3. This simulation concludes the case 
study for the energy space E defined in Section 2.2. 

Now, if we want to use the generalized energy operators (and Theorem 2) with p 
in Z+ and p > 0, then one can look at solutions in the energy space Ep (e.g. Section 
2.2). It is the solutions of the kind ([[/]^]i'’*)"' or {[[f]^]i’^)^. Thus, the definition 
of the subspaces and N 2 can be extend with p in Z+ such as 

Nyp = { 5 eS"(K^), Vte[0,T], ro e [ri,r 2 ]| 5 (ro,t) = at*([[/]P])^’*)”(ro,t), 

/ € Sp (M^), n G Z"*" — {0,1}} 

N 2 ,p = {g& S~(]R^),Vr e [ri,r 2 ], to G [0,T]|g(r,to) = 9);([[/]*’]+’’')"(r,to), 

/ G Sp (M^), n G Z"*" — {0,1}} 

With this definition, let us do a similar numerical estimation of the solutions and 
finding ii and 12 for the case p = 1. The bottom figures of Figure 1 and 2 are 
the results. We can see that there is a factor of 10“'* and 10“ 11 in Scenario 1 
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Figure 2. Results from the estimation of the energy of 
and 9®<o))- The Figures correspond to 

Scenario 2. 



>. 

o> 


LU 

■o 


£ 



i-th Derivatives ° 




and 2 respectively between the estimated energy for the operators corresponding 
to the case p = 0 and p = 1. Thus, one can conclude that ii and are equal to 
0. Recalling the decomposition Theorem 1 (and equation ( |3.3| )), the energy space 
Eo is reduced to {0}. 


4.3. Phenomenon of Duplication of Waves in a Closed Cavity. This section 
reflects on the idea of propagating electromagnetic waves in a closed cavity [2]. Here, 
the idea is not to revise the whole theory of electromagnetism, but rather discussing 
at a theoretical level the implication of our proposed model in this example and 
defining the duplication of waves. 

Let us consider the form of solutions which propagates in a closed cavity. According 
to m, one possible solution is the evanescent wave described in equation ( |4.4[ ). One 
can notice that / and £{f) are analytic by definition for the variable t and r. With 
Proposition 1, we can assume that / is finite energy and in s“(IR^) with a wise 
choice on the parameters A, fci, k 2 and w. One can estimate the difference of energy 
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in time over dt inside the cavity at a specific location rg (rg in [ri,r 2 ]) such as 


[ {f{ro,u)fdu < oo 
Jo 


£ifiro,T)) 

£{f{ro,T + dt)) = £{f{ro,T)) + Y,d^if{ro,T))^<oo 

k^O 

£{f{ro,T + dt)) = £:(/(rg,r)) + ^a^=(f(rg,T))^<oo 

oo 

£{f{ro,T + dt)) = £{f{ro,T)) + f{ro,T)dt + J2dt\'i’V{f)iro,T) 

+vI/-‘(/)(rg,T))(^<cx. 

£{f{ro,T + dt)) ~ £{f{ro,T)) + firo,T)dt (4.8) 

Note that by definition dt)"’*(/) (rg, T) = 0. Here the symbol means that 

3eGM+, e«l, Vk G Z+ , k > 0\ \dt\^t’\f)iro,T))\ < e\f{ro,T)\ (4.9) 

Now, let us do a hypothesis that £(/(rg,T + dt)) increases significantly over dt 
modifying the approximation in ( |4.9[ ) 

3eGK+, e«l, Vfc G Z+, k > 1\ \dt^'i>t'\f)iro,T)\ < e\'i't'\f){ro,T)\ 

(4.10) 

and then, 


£ifiro,T + dt)) ~ £{firo,T)) + firo,T)dt + ^t’{f)iro,T) 


dt^ 


(4.11) 

Now, using the Theorem 1 and the model based on the energy space in the previous 
sections, one can say that the subspace N( defined in Section 4.1 is not reduced 
to {0}. Thus, the solutions from this subspace has to be taken into account. The 
duplication of wave can be formulated as an approximation for taking into account 
additional solutions produced by the wave equation. 


5. Conclusions 

The core of this work is to define the notion of multiplicity of the solutions of 
a linear PDE using the model associated with energy spaces and the (generalized) 
energy operators. In this way, it contradicts the classical way of solving a nomi¬ 
nated PDE with boundaries conditions, but it rather focuses on additional solutions 
from these energy spaces. The multiplicity is defined through Theorem 2 and the 
Corollary 1. The work shows how the energy operators (and generalized energy 
operators) can determine which energy subspace is reduced to {0}. 

The theory is then applied to the evanescent waves, a special type of solutions of 
the wave equation. The last part with the closed cavity show a possible real world 
applications. In this case, the duplication of waves is when additional solutions 
should be taken into account due to the level of energy increasing in the cavity. 
However, this work remains at a theoretical level and more work with simulations 
are required to fully understand the concept of duplication. 
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Furthermore, it was explained in Section |3.4| that the linearity of the PDFs does 
not play any role in Theorem 2 and the Corollary 1. As a matter of fact, the 
model with the energy subspace is applied directly to the form of the solutions of 
the nominated PDF (i.e. evanescent waves solution of the wave equation) with the 
conditions that the solutions are finite energy and in the Schwartz space S“(IR^). 
Thus, there is no restriction theoretically speaking to use this model with non¬ 
linear PDFs. Our next interest is to apply our model to the type of solutions of the 
Korteweg de Vries equation called solitons which filled in the properties required to 
apply this model (e.g, finite energy function, solution in the Schwartz space S~(IR^) 
)• 

To conclude, the multiplicity of the solutions applied to the wave equation is only 
the first step before wondering if the receiving waves in a transmitter-receiver sys¬ 
tem, can be a sum of all those solutions. In other words, one can ask if the signal 
is generated by receiving not one type of solution/wave, but the additional solu¬ 
tions/waves coming from other energy subspaces. 
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Appendix A. Discussion about n in Z 

Here is the discussion written in [15] (Section 4) in the case of the generalized 
energy operators about extending Theorem 1 and Lemma 1 with n in Z. The 
extension of Theorem 0 and Lemma 0 can be found in [14] (Section 3, p. 74). 
Let us recall: 

Discussion u < — 1: In this case, one can define: 


v/es- 


VteR, peZ+, {[[f{t)]P]tP ^0, VnGZ+,n>l, 


mmtr 

(A.l) 


This set of functions can also be described as: / in Sp 
for p in Z+. Note that one could also chose to have / in Sp 


and / not in Ker[[[f{t)]P]^) 


However, this is 


more restrictive than the set defined in (A.l). Using an intermediary function, h 
such as h = ; the problem of decomposing i9^([[/(t)]^]i’) " {k in Z+ — {0}) 

is equivalent to resolving which has been demonstrated in the Lemma 1 and 

Theorem 1. 


Discussion n = loru = — 1: As already underlined in [14], one can use a general 
formula for / in the set defined in equation 


dHiimW 

k = l, 5*([[/(f)f]+) 
k = 2, 


‘ v([[/(t)M+)V 

{[[fmty'mmW + {[imntfmfmty" 
+mmtYdy[[fmty^ + {[[mryy'dunfmtf 

(A.2) 


The example for k = {1,2} in Equation (A.2) shows that dt {[[f{t)]^]t) can be de¬ 


composed into a product of successive derivatives of ([[/(O]^]!")^ ([[/(t)]^]i~) 

Those derivatives can be decomposed into a sum of generalized energy operators 
based on the Lemma 1 and Theorem 1 plus the previous discussion (for the case 
n < —1 ). 








MULTIPLICITY OF SOLUTIONS FOR LPDES 


17 


Now for the case n = — 1, it is easy to see that : 


dHilfimty' 




(A.3) 


With the discussion for the case n = 1, we can conclude that 9(^([[/(t)]^]i') ^ 
can be decomposed into a product of successive derivatives of ([[/(t)]^]^)^ and 

{{{fmty'- 
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